Current and vortex statistics in microwave billiards 
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Using the one-to-one correspondence between the Poynting vector in a microwave billiard and the 
probability current density in the corresponding quantum system probability densities and currents 
were studied in a microwave billiard with a ferrite insert as well as in an open billiard. Distribution 
functions were obtained for probability densities, currents and vorticities. In addition the vortex 
pair correlation function could be extracted. For all studied quantities a complete agreement with 
the predictions from the random-superposition-of-plane-waves approach was obtained. 
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PACS numbers: 05.60.Gg, 05.45.Mt, 73.23.Ad 

I. INTRODUCTION 

A particularly successful approach to describe the sta- 
tistical properties of the eigenfunctions of chaotic bil- 
liards, dating back to Berry tjj|, assumes that at any point 
not too close to the boundary the wave function can be 
described by a random superposition of plane waves, 



tk n r 



(1) 



where the modulus k = \k n \ of the incoming waves is 
fixed, but directions k n /k and amplitudes a n are con- 
sidered as random. In billiards with time-reversal sym- 
metry there is in addition the restriction that the wave 
function has to be real. This ansatz cannot be strictly 
true. It completely ignores the boundary conditions at 
the billiard walls, but this is of no importance as long as 
the wave length is small compared to the billiard size. 

As a consequence one expects a Gaussian distribu- 
tion for the wave function amplitudes tp, or, equiva- 
lently, a Porter-Thomas distribution for their squares 
\tp\ 2 . For the spatial auto-correlation function c(r) = 
(i/j*(f + r)if)(f)) f a Bessel function is found. These pre- 
dictions were first verified by McDonald and Kaufman 
in their disseminating papers on stadium wave functions 

0]. It is impossible to mention all works which have 
been published hitherto on this subject. As another con- 
sequence of ansatz (P the gradient of the wave function is 
Gaussian distributed, too, and uncorrelated to the wave 
function. This has been used to calculate the distribution 
of eigenvalue velocities and the velocity auto-correlation 
function for the case of a local parameter variation Q . 

The approach is not restricted to quantum-mechanical 
systems, or systems where there is a one-to-one cor- 
respondence to quantum mechanics such as quasi-two- 
dimensional microwave billiards. Porter-Thomas distri- 
butions are found as well in the squared amplitudes 
of vibrating plates JBj. In an experiment on a three- 
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dimensional Sinai microwave billiard, having no quan- 
tum mechanical equivalent, the electromagnetic field dis- 
tributions as well as their spatial correlations could be 
explained assuming a random superposition of plane elec- 
tromagnetic fields || 0]. In a very recent experiment on 
light propagation through distorted wave guides, finally, 
the observed transversal field patterns could be described 
again by ansatz ([j]) || . All these examples show that the 
random-superposition-of-plane-waves approach is not of 
a quantum-mechanical origin, but holds for all types of 
waves. 

If the billiard is opened, or if time-reversal symmetry 
is broken, the wave function acquires an imaginary part, 

with the consequence that in dependence of the relative 
fractions of real and imaginary parts the distribution of 
\i(j\ 2 changes from Porter-Thomas to single exponential 
behavior. An explicit formula describing the distribution 
in the transition regime has been given by different au- 
thors H [Uj [ll], [l^] . The same function has been derived 
by Seba et al. for the distribution of scattering matrix el- 
ements in a partially opened microwave billiard (T3J. Wu 
et al. studied amplitude distributions and spatial auto- 
correlation functions in a microwave billiard with one 
ferrite-coated wall to break time-reversal symmetry, and 
found quantitative agreement with the results expected 
from the random-superposition-of-plane-waves approach 
ifTijl . In a recent paper by Ishio et al. also deviations 
of this formula due to scars and in regular systems are 
studied |||. 

More recently the interest focussed on the current 
statistics in open systems. Saichev et al. calculated the 
distribution of currents p2| . The properties of current 
vortices have been studied by Berry and Dennis, who 
gave analytic expressions for different types of vortex spa- 
tial auto-correlation functions Independently, such 
auto-correlation functions as well as the distribution of 
nearest distances between nodal points have been studied 
in Ref. @. 

Experimental results on current statistics and the dis- 
tribution of nodal points have not been available as yet. 
The method of choice to study such questions are mi- 
crowave techniques. In a previous work of our group we 



2 




FIG. 1: Map of I'i/'l 2 (top) and of the current (bottom) in the ferrite billiard at 5.19 GHz (a) and 6.41 GHz (b), and for the open 
billiard at 5.77 GHz (c). The intensity has been converted into a gray scale with the intensity increasing from white to black. 
The lengths of the arrows in the current map correspond to the magnitude of the Poynting vector. For the field distribution 
in panel (a) there is a strong flow into the probing antenna at points of high intensity. In the two other examples the flow into 
the probing antenna is negligible. 



could demonstrate that currents can be easily obtained 
in an open microwave billiard [ fl8| . In the present paper 
we show results on current distributions in two types of 
billiards, verifying a number of predictions given in the 
above mentioned papers for the first time also in experi- 
ment. 



II. EXPERIMENT 

The basic principles of the experiment are described 
in detail in Ref. Therefore we concentrate on the 

aspects being of relevance to the present study. Two dif- 
ferent billiards have been used, both resonators have a 
size of about 25 cm and a height of 8 mm. One of them 
is an open system of a rectangular shape with rounded 
corners and two openings on opposite sides containing 
entrance and exit antenna (see Fig. Et). A third movable 
antenna was used to map the field distribution in the 
resonator on a square grid of period 5 mm. The same 
system was used for a quantum-dot analog study previ- 
ously ]|(|. The two fixed antennas have a metallic core 
of diameter 1 mm, and a teflon coating for stronger cou- 
pling; the movable antenna was a thin wire of diameter 
0.2 mm. The lenghts of all antennas were about 6 mm. 



The quantity directly accessible in the experiment is the 
scattering matrix ]2l"| ] . From a reflection measurement as 
a function of the antenna position a mapping of the mod- 
ulus of the wave function can be obtained; to get the sign 
as well, the transmission between two antennas is needed. 
By the presence of the antennas the resonances are some- 
what broadened and shifted, but in chaotic systems this 
does not change the universal distributions. (For details 
see Ref. @.) 

The second system studied is a Robnik, or limacon bil- 
liard |22|, |2^] . It can be obtained by a complex mapping 
of the unit circle in the complex plain by means of the 
function w = z + \z 2 . The deformation parameter used 
in the experiment was A = 0.4. For this parameter the 
billiard is believed to be completely chaotic. (There are 
recent theoretical results indicating that this cannot be 
taken for sure ^ij. But if there are still stable islands 
in the phase space, their corresponding volume will be 
very small.) A ferrite ring is introduced to break time- 
reversal symmetry (shape and position are included in 
Figs. [l]a and b). Ferrites have been used already repeat- 
edly for this purpose @ ||, H' ■ ^ detailed descrip- 
tion of the function principle of the ferrites will be given 
in a forthcoming publication 28 . During reflection at 



the ferrite the microwaves experience a phase shift with 
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the consequence that there will be currents through the 
billiard though it is completely closed. These currents 
are a complete analogue to the persistent currents ob- 
served in mesoscopic structures P9| , for a recent review 
see Ref. [3C|. They will be the subject of a separate pub- 
lication 1 31]. In addition there will be currents due to 
the fact that the ferrite introduces considerable absorp- 
tion into the system. In the present context it is not of 
relevance whether the currents are due to a break of time 
reversal symmetry or due to absorption. 

It was demonstrated in Ref. |Q that in quasi-two- 
dimensional billiards the Poynting vector S = -f-E x H 
can be written as 



8?rfc 



Im [£*(r)V£(r)] 



showing that there is a one-to-one correspondence with 
the current density 

j = — Im [?/>*(r)V?/>(Y)] 

in the corresponding quantum-mechanical system. A 
measurement of the electric field in the resonator includ- 
ing the phase thus immediately yields the Poynting vec- 
tor and by means of the mentioned analogy the current 
density. 

There is one problem with the experimental determi- 
nation of field and current distributions. The probe an- 
tenna moving through the billiard unavoidably gives rise 
to a leakage current spoiling the statistical properties of 
the current distribution pq |. The influence of the probe 
antenna is small as long as there is a strong flow through 
the system. In this case the unavoidable leakage current 
into the probe antenna is negligible. There are situations, 
however, where there is no or only little flow, e.g. for the 
open dot, if the total transmission is close to zero, or 
for the ferrite billiard, if there are strong standing waves 
present, as in Fig. |l|a. In such cases the leakage current is 
no longer negligible, and may even become dominating. 
There is only one way to avoid this problem: the frequen- 
cies have to be chosen such that the overall amplitudes 
are moderate. As an example Fig. [l]a shows a mapping of 
\ip\ 2 and of the current distribution in the ferrite billiard 
at a frequency where there is a strong flow into the mea- 
suring antenna at points of high intensity. In Fig. 0b, on 
the other hand, the frequency has been chosen such that 
there is no noticeable current into the probe antenna. 
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FIG. 2: Distributions of p = |i/>| 2 for the ferrite billiard at 
5.19 GHz (a), 6.41 GHz (b), and for the open billiard at 
5.77 GHz (c), corresponding to Fig. [j]. The dashed lines have 
been calculated from Eq. 



III. INTENSITY DISTRIBUTIONS 

In a chaotic billiard real and imaginary part of the 
wave function are uncorrelated within the random-super- 
position-of-plane-waves approach, (ipRipi) — 0, where the 
average is taken over the billiard area. In a microwave 
measurement, there may be an additional rotation of the 
phase, caused by antennas, cables etc., i.e. the wave 



function registered by the apparatus is 

i> = 4>r + ii>i = e ia {ip R + iipi). 

This has the consequence that real and imaginary part 
become correlated, (tpRtpi) ^ 0. To begin with we have 
to remove this phase rotation. The procedure is de- 
scribed, e.g., in Refs. (l^, |Hj |. 

The distribution of intensities p = \tp\ 2 can be calcu- 
lated with help of the Berry ansatz described above. If 
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FIG. 3: Distributions of j x (a), j y (b), and j, the latter one both on a linear (c) and a logarithmic scale (d), for the ferrite 
billiard at 6.41 GHz, corresponding to Fig. ^b. The dashed lines are the theoretical expectations from Eqs. (|^) and (Q), 
respectively, where a — yj (jx,y), see Eq. (|^). 



there is a complex mixing of independent Gaussian fields 
ijjR and ipi, we have in the transition regime |12l n5[ 



P(p) = (j, exp(-^ 2 p)^o( ) "V / Ai 2 - lp) 



(2) 



where 



and 



'MI 



In the limit e — > 0, Eq. (g) describes the Porter-Thomas 
distribution found for systems where the wavefunction is 
real. For e — > 1, a single exponential behavior is observed 
which holds for systems where real and imaginary part 
are of the same strength. 

Figures ^a and b show the corresponding intensity dis- 
tributions for the ferrite billiard. The dashed lines are 
calculated from Eq. (0). The parameter e has not been 
fitted but has been taken directly from the experiment by 
averaging ijj R and ipj over all pixels at a given frequency. 

One observes a quantitative agreement with theory for 
the case that the leakage current is negligible (Fig. ||b). 
For the scarred wavefunction of Fig. @a, on the other 
hand, the experimentally found distribution of |?/>| 2 is 



completely at odds with theory. Such discrepancies are 
not new. Already in the disseminating paper by McDon- 
ald and Kaufman || Gaussian distributions for if> were 
observed for chaotic wave functions exclusively. It is ob- 
vious that the random-superposition-of-plane-waves ap- 
proach cannot work for bouncing-ball and scarred wave 
functions such as the one shown in Fig. l|a. It is an easy 
matter to show that sharp drops in the ip\ 2 distribution 
as the one in Fig. ||a are generic, e.g., for all wave func- 
tions of the rectangle 

Figure |^c shows a corresponding example for the open 
billiard at a frequency where there is no leakage current. 
Again one finds a complete correspondence with theory. 



IV. CURRENT DISTRIBUTIONS 

For the currents shown in the bottom panel of Fig. |l| 
the distribution function of 

j = Im [ip*Vip] = ^PrV^i - tpiVtpR 

has to be calculated. A similar problem occurred in our 
work on global and local level dynamics where the dis- 
tributions function of ipVip was needed B. With the 
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FIG. 4: Same as Fig. ^[ but for the open billiard at 5.77 GHz, corresponding to Fig. [j]:. 



random-superposition-of-plane-waves approach the cal- 
culation of the averages is straightforward and yields |12] 



PU*,v) = 



1 



201,,) 



: exp 



(jx.y) 



\jx,y I 



(3) 



for the distribution of the current components j x and j y , 
and 



P(j) = %-K { 



2j 



0' 2 > "VVO' 2 ), 

for the distribution of the modulus j = 
the parameter 



(4) 



Jx + i« where 



(5) 



can again be taken directly from the experiment. 

Figures |3| and ^ show the distributions for j x , j y , and 
j for the same patterns depicted in Figs. |l]b and c, re- 
spectively. Once more theory is in complete accordance 
with the experiment if there are no leakage currents, but 
is not able to describe the experiment otherwise. For 
comparison, Fig. ^ shows the current distribution for the 
scarred wave function in Fig. |l|a. Not unexpectedly there 




FIG. 5: Distribution of j on a logarithmic scale as in Fig. gd, 
but for the scarred wave function at 5.19 GHz, shown in 
Fig. 0a. 



are again significant deviations from the universal behav- 
ior. The excess at large current values in particular is a 
consequence of the strong leakage currents observed at 
positions of high |?/>| 2 values (see Fig. |l]a). 

In particular we observe identical distributions for j x 
and j y in the ferrite billiard. This is no longer the case 
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for the open billiard, where the maximum of the j x dis- 
tribution is shifted significantly to negative values caused 
by the transport from the right to the left through the 
billiard. 



V. THE VORTEX SPATIAL 
AUTO-CORRELATION FUNCTION 

A very useful quantity to characterize the vortex struc- 
ture of a vector field is the vorticity. Up to a factor 1/2 
it is identical with the rotation of the current density, 



1.00 



iv = 



dx dy dx dy 



Vx ] 



The calculation of the corresponding distribution func- 
tion P{uj) follows exactly the same lines as for the current 
distribution and yields 



i 



exp 




where 



<^ 2 > = Wu&Htf)- 



(6) 



(7) 



can be taken directly from the experiment. 

Figure ^ shows the vorticity distribution for the exam- 
ple shown in Fig. |l|b for the ferrite billiard. The dashed 
line has been calculated from Eq. (^). Again the param- 
eter A was not fitted but taken from the experiment. 

For the wave function to be zero both real and imag- 
inary part have to be vanish. As a consequence there 
are no longer nodal lines, but only nodal points. Each 
nodal point corresponds to a vortex in the correspond- 
ing flow pattern. Since the distance between neighboring 
node lines is of the order of half of a wave length both 
for real and imaginary part, the mean spacing between 
neighboring nodal points is of this order of magnitude as 
well. 

Using the random-superposition-of-plane-waves ap- 
proach, Berry and Dennis jl6| calculated two types of 
vortex spatial auto-correlation functions, one of them 
called pair correlation function and defined by 

9{r) = g {5[iljR{f + r)]S[ipi(f + r)]6[ip R (f)]S[^i(f)] 
x | w (f + r)|Kf)|) f . (8) 

The normalization go is determined such that g(r) — > 1 
for r — > oo. The other correlation function discussed by 
the authors, the charge auto-correlation function gq{r), 
differs from expression (^) only by the fact that the sign 
of the vorticity is taken into account, 

9Q<X) = 9o(S[il) R (r + r)]5[<ipi(f+r)]5[il) R (r)]5[<ilJr(r)] 
x w(f + r) u>(f)) f . 

Since for gq{r) pairs of vortices with different senses of 
rotation enter with a negative sign, we have gQ{r) — > 
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FIG. 6: Distribution of the vorticity for the ferrite billiard at 
6.41 GHz corresponding to Fig. lib. The dashed line corre- 
sponds to Eq. (^|), where A = (w 2 ), see Eq. (j?j). 




FIG. 7: Vortex pair correlation function g(x) for the ferrite 
billiard, see Eq. (^), where x — kr and g(x) = g(r). The solid 
line has been calculated from an integral expression given by 
Berry and Dennis [M. 



for r — > oo. For the explicit expressions of g{r) and gq{r), 
which are quite complicated, the reader is referred to the 
original work. From the spatial auto-correlation function 
the distribution of nearest distances between vortices can 
be calculated which has been studied by Saichev et al. 

0- 

The experimental study of the vortices is more difficult 
than that of the different types of distribution functions 
discussed above, since it involves the determination of 
the zeros of real and imaginary part of the wave func- 
tion. Because of the applied grid period of 5 mm, the 
precision in determining distances is only moderate, and, 
more problematically, it is impossible to resolve vortices 
lying very close together which is very often the case for 
vortices with different sign. In this situation a reliable 
determination of the distance distribution is not possi- 
ble. The determination of the charge correlation func- 
tion <7q, too, did not work because of the cancellation of 
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positive and negative terms, leading to an intolerable in- 
crease of the noise level. But the pair correlation function 
g could be determined. Since the wave number k enters 
the pair correlation function as a scaling factor only, we 
may write g(r) = g(x), where x = kr. This allows to 
improve the statistics by superimposing the results from 
different frequencies. Figure ^ shows the resulting pair 
correlation function g{x) obtained from current patterns 
at 43 different frequencies in the range from 4 to 8 GHz, 
by extracting the found vortices by hand. The frequency 
regimes showing a flow into the probe antenna were care- 
fully avoided. Though the statistics is only moderate, 
the oscillations predicted by theory are clearly observ- 
able. The hole observed in the experimental histogram 
at small distances reflects the above mentioned difficulty 



to resolve closely neighboring vortices. 

We can thus conclude that the experimentally ob- 
tained distributions of wave function amplitudes, cur- 
rents, and vortices are in quantitative agreement with the 
predictions of the random-superposition-of-plane-waves 
approach, if the wavefunction is chaotic. 



Acknowledgments 

Discussions with K.-F. Berggren, Linkoping, and with 
M. Vranicar, M. Robnik, Maribor, are gratefully acknowl- 
edged. The experiments were supported by the Deutsche 
Forschungsgemeinschaft . 



[15] 



M.V. Berry, J. Phys. A 10, 2083 (1977). [16 
S.W. McDonald and A.N. Kaufman, Phys. Rev. Lett. 42, 
1189 (1979). [17 
S.W. McDonald and A.N. Kaufman, Phys. Rev. A 37, 
3067 (1988). [18 
M. Barth, U. Kuhl, and H.-J. Stockmann, Phys. Rev. 
Lett. 82, 2026 (1999). [19 
K. Schaadt, Master's thesis, Niels Bohr Institute, Copen- 
hagen, 1997. [20 
U. Dorr, H.-J. Stockmann, M. Barth, and U. Kuhl, Phys. 
Rev. Lett. 80, 1030 (1998). [21 
B. Eckhardt, U. Dorr, U. Kuhl, and H.-J. Stockmann, 
Europhys. Lett. 46, 134 (1999). [22 
V. Doya, O. Legrand, F. Mortessagne, and C. Miniatura, [23 
Phys. Rev. Lett. 88, 014102 (2002). [24 
K. Zyczkowski and G. Lenz, Z. Phys. B 82, 299 (1991). [25 
E. Kanzieper and V. Freilikher, Phys. Rev. B 54, 8737 
(1996). [26 
R. Pnini and B. Shapiro, Phys. Rev. E 54, R1032 (1996). [27 
A.I. Saichev, H. Ishio, A.F. Sadreev, and K.-F. Berggren, 
J. Phys. A 35, L87 (2002). [28 
P. Seba et al., Phys. Rev. E 56, 2680 (1997). [29 
D.-H. Wu, J.S.A. Bridgewater, A. Gokirmak, and S.M. 
Anlage, Phys. Rev. Lett. 81, 2890 (1998). [30 
H. Ishio, A.I. Saichev, A.F. Sadreev, and K.-F. Berggren, [31 
Phys. Rev. E 64, 056208 (2001). 



M.V. Berry and MR. Dennis, Proc. R. Soc. Lond. A 456, 
2059 (2000). 

A.I. Saichev, K.-F. Berggren, and A.F. Sadreev, Phys. 
Rev. E 64, 036222 (2001). 

P. Seba, U. Kuhl, M. Barth, and H.-J. Stockmann, J. 
Phys. A 32, 8225 (1999). 

U. Kuhl, E. Persson, M. Barth, and H.-J. Stockmann, 
Eur. Phys. J. B 17, 253 (2000). 

Y.-H. Kim, M. Barth, H.-J. Stockmann, and J. P. Bird, 
Phys. Rev. B 65, 165317 (2002). 

J. Stein, H.-J. Stockmann, and U. Stoffregen, Phys. Rev. 
Lett. 75, 53 (1995). 

M. Robnik, J. Phys. A 16, 3971 (1983). 

M. Robnik, J. Phys. A 17, 1049 (1984). 

HR. Dullin and A. Backer, Nonlinearity 14, 1673 (2001). 

P. So, S.M. Anlage, E. Ott, and R.N. Oerter, Phys. Rev. 

Lett. 74, 2662 (1995). 

U. Stoffregen et al, Phys. Rev. Lett. 74, 2666 (1995). 
H. Schanze, E.R.P. Alves, C.H. Lewenkopf, and H.-J. 
Stockmann, Phys. Rev. E 64, 065201(R) (2001). 
H. Schanze et al, to be published. 

L.P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, 

Phys. Rev. Lett. 64, 2074 (1990). 

P. Mohanty, Ann. Phys. (Leipzig) 8, 549 (1999). 

M. Vranicar et al., J. Phys. A (in press) 



